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Abstract 

We investigate what one can hope to learn about the parameters that describe the neutrino 
fluxes emitted by the explosion of a galactic supernova using the observations of a megaton-size 
Water-Cherenkov detector. We calculate the allowed regions that can be obtained by fitting these 
parameters to a simulated observation of events by such a detector. All four available detection 
channels (inverse beta decay, charge and neutral current on oxygen and elastic scattering on elec- 
trons) are included in the fit and we use a ten dimensional parameters space. Nine parameters 
describe the initial neutrino fluxes and are referred to as the supernova parameters. Furthermore, 
we include the dependence on the Chooz mixing angle ^13, which controls the matter effects that 
the neutrino undergoes in the outer-parts of the supernova. If we do not make any extra assump- 
tion on these parameters, we show that one can hope to determine ^13 quite well whereas, except 
for the parameters describing the Ve flux, most of the supernova parameters are rather difficult to 
constrain, even if the four detection channels could be completely separated. 
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I. INTRODUCTION 



The phenomenon of a core-collapse Supernova (SN) explosion embraces vastly different 
areas of physics and involves all the known fundamental interactions The core-collapse 
itself is driven by gravitational effects; the thermodynamics is controlled by electromagnetic 
and strong forces and naturally the weak force plays a major role in the energy loss through 
emission of neutrinos. Evidently, the study of supernovas is a promising playground for 
testing new physics and to learn about particle properties. 

Among the particles one hopes to learn more about are the neutrinos. Indeed the major 
part of the enormous energy liberated in a core-collapse (type-II) supernova is emitted in 
the form of neutrinos. Naturally, one might also attempt to use our knowledge of neutrino 
properties to extract information about the supernova physics from observations of the 
neutrino fluxes. In particular, as the neutrinos are emitted from the interior of the star, 
they offer a way to probe the physics of the core-collapse. But independent of whether 
one attempts to learn about supernova physics or about neutrino properties from a SN 
observation, ultimately, when comparing data with theory one needs to take into account 
the uncertainties on the supernova parameters as well as on the neutrino parameters. 

Many efforts have been made in order to understand the complex physics involved and to 
predict the emitted neutrino fluxes . When discussing supernova neutrino there are only 
three distinguishable flavors, as the fi and r neutrinos, as well as their antineutrinos, have 
identical properties. We will denote the /i and r neutrinos and antineutrinos by u^. Thus 
along with the electron neutrino (z/e) and the electron antineutrino (z/g) constitute the three 
supernova neutrino species. The neutrinos from the core-collapse supernova are expected 
to be almost thermal, since they are in fact trapped inside what is known as the neutrino- 
sphere. Besides the very early universe a supernova is probably the only place neutrinos 
thermalize. However, small deviations from the thermal spectrum are expected and how 
narrow or wide the energy spectra will be is in general described by a pinching parameter. 
Therefore, it is a good approximation to parametrize each initial neutrino spectrum by three 
parameters; the average neutrino energy, the total emitted energy (or integrated luminosity) 
and the pinching parameter. 

During the last few years evidences of neutrino masses and mixing have been gathered and 
today most of the neutrino oscillation parameters have been determined with good precision, 
from the observation of solar and atmospheric neutrinos as well as other terrestrial neutrino 
experiments 

aaaaa. 

Some important parameters remain unknown. The absolute 
neutrino mass scale, mo, is only bounded from above by cosmological data mo ~ 0.2 — 0.7 
eV jol. The same is true for the mixing ang le 6*13, sin2 6'i3 < 0.04 0. Furthermore, the 
pattern of neutrino masses is not established yet: we do not know if nature has chosen the 
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normal {m^ > m2 > mi) or the inverted (m2 > rrii > m^) mass hierarchy, where mi (m^) is 
the mass of the neutrino state most (least) populated by the z/g component. Finally, we do 
not know anything about the existence of CP violating phases in the neutrino sector *. 

In particular, for supernova neutrinos the neutrino mass hierarchy and the Chooz angle, 
^13, are very important, as the so-called Mikeyev-Smirnov-Wolfenstein (MSW) H resonance 



strongly depends on them 12|, [iSj. Recently, there has been several investigations of the 
possible impact of the neutrino-neutrino interactions in the dense neutrino region inside the 
neutrino-sphere in the supernova [3]. These self interactions may cause collective effects 
and thus influence the neutrino survival probability as a function of energy. Under realistic 
supernova density profiles the collective effects are negligible for the normal mass hierarchy. 
In the case of the inverted hierarchy it seems that a swapping of the Ue and energy 
spectra above a certain critical energy occurs, and at the same time a swapping of the 
antineutrino spectra occurs. In the latter case there are indications that the spectra might 
be smeared out and therefore does not exhibit a sharp interchange of the spectra at a certain 
energy. In principle, these swapped spectra should then be used as input when calculating 
the MSW effects in the outer parts of the supernova. However, since the impact of the 
neutrino-neutrino self interactions is still being debated, we will in this paper only consider 
the normal mass hierarchy and we can thus neglect the collective neutrino-neutrino effects. 

One can imagine two very different sources of supernova neutrino signals. The diffuse 
supernova neutrino background (DSNB), arising from all past supernova explosions; and 
the lucky event of a galactic supernova. In the latter case one expect about 10^ — 10^ 
neutrino events in a megaton scale detector in a 10 seconds time interval. Therefore, this 
is practically free of background. On the other hand, background is the major issue for the 
DSNB detection and one expect only to be able to extract a signal in a rather small energy 
window and most likely only for the dominant detection channel. Therefore, depending on 
the detector, there will be sensitivity only to the electron antineutrino flux (Water- Cherenkov 
and Scintillator) or only to the electron neutrino flux (Liquid Argon). On the contrary, for 
galactic SN observations one expect to have several detection channels available, each with 
sensitivity to different neutrino flavor compositions - a feature which is crucial for the pinning 
down of the neutrino parameters. As is well-known, a galactic supernova explosion is a rare 
event. About two per century is the best we can hope for. The optical observation is actually 
likely to be obscured by dust etc., and in this light it seems even more crucial to get a better 
understanding of the SN neutrino fluxes. 

Here we will focus on the investigation of what can be learned from an observation of the 



CP-violation effects are ignored in our analysis. It was explicitly shown in llj that there is no net CP 
effect due to standard neutrino oscillations in a SN core-collapse. 
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neutrino burst from a galactic supernova. We will furthermore restrict ourself to what can 
be measured by a Water- Cherenkov detector. We aim at figuring out what can be learned 
about supernova physics as well as what can be learned about the neutrino parameters simul- 
taneously. We thus continue previous studies 15|, ll6|, ll7| , but we will increase the parameter 
space, taking into account all important parameters in what optimistically could be referred 
to as a full parameter space fit. In the previous studies only the dominant detection channel 
has been considered or some of the parameters have been fixed. In most earlier works the 
pinching parameters were fixed, but in Ref. 17j it was shown that, when considering only 
the dominant inverse beta decay detection channel, these constitute important uncertainty 
factors when attempting to extract information about the supernova parameters. 

In fact, the authors of Ref. 17| have discovered a degeneracy, which only appears when 
including the pinching parameters, between the mean energy (Ex), the luminosity and 
the pinching parameter f]x of flux. By degeneracy we mean that, for a given allowed point 
in {{Ex), Lx, (3x)i i.e. a point where the data are explained (one can always take the input 
values themself ) , it is possible to construct a different allowed point which has a larger value 
of {E^) by increasing (3x and decreasing the luminosity L^- As this degeneracy covers almost 
the complete parameter space, the determination of the correct value of {E^) (as well as (S^ 
and Lx) becomes very hard. Indeed the allowed region for {Ex) was shown to include the 
whole range from 15 MeV to 30 MeV as expected from supernova simulations 17|]. This 
degeneracy occurs since shifting the parameters in the mentioned fashion, maintains the 
energy-spectrum at the Earth almost identical, with only small differences in the low energy 
spectrum. Due to the threshold (which is about 5 MeV) on the inverse-beta decay channel, 
the Water- Cherenkov detector is not sensitive to the spectrum at low energies. Clearly 
the Ve parameters have also to be adjusted correspondingly, but as the allowed regions are 
still rather small, we will not speak of a degeneracy in these variables. However, only the 
inverse beta decay channel has been considered in Ref. Q and .fsuot dear to what extend 
the inclusion of the other detection channels might break down this degeneracy. Indeed in 
Ref. [isl several degeneracies in the only-inverse-beta-decay case were found to be broken 
down when considering all four detection channels. This is one of the motivations for the 
present work. 

In our full parameter space fit, we use a total of ten fitting parameters. Nine parameters 
describe the initial supernova neutrino flux (z/g, Ve and luminosities, mean energies and 
pinching parameters), as stated above three parameters are needed for each of the three 
neutrino species. Furthermore, we also freely vary the most important neutrino parameter, 
the Chooz angle ^13. However, as explained above we fix another important neutrino param- 
eter, the neutrino mass hierarchy, as normal. The terminology full parameter space should 
therefore, as usual, be taken with a grain of salt. There might be a number of other pa- 
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rameters, such as deviation of the supernova density profile from the assumed p ^ form [1 
effects from sterile neutrinos 



20|] or from new interactions [ll|, deviation from the assumed 
spectral forms and yet undiscovered effects that may influence the flux of neutrinos from a 
supernova. 

II. THE ANALYSIS PROCEDURE AND THE PARAMETER SPACE COVER- 
AGE 

In this section we describe the details of our analysis procedure and the chosen parameter 
space. As mentioned in the introduction we will investigate how well the supernova and 
neutrino parameters can be determined by the observation of neutrinos from a galactic SN 
by a future megaton scale Water-Cherenkov detector. We will take into account all major 
parameters (with exception of the neutrino mass hierarchy, here assumed to be normal) that 
can influence the neutrino flux, thus expanding the parameter space as compared to earlier 
works. We vary a total of 10 parameters to be fltted by data: the Ve-, and luminosities 
(Le, Le and Lx), mean energies {{Ee), {Eg) and {E^)), pinching factors (/3e, Pe and Px), 9 
SN parameters, and a single neutrino quantity, sin^6'i3. 

We will assume that the initial supernova neutrino fluxes emitted at the respective 



neutrino-spheres can be parametrized with the spectrum as suggested in Ref. 21|. Thus 



for each neutrino specie i = e,e, x, we assume the energy spectrum to be of the form 

/ E Y'-' , ^ E , 



exp(-/3j- 



where (3i is the pinching parameter, E the neutrino energy and (Ei) the Ui mean energy. 
For Pi > 3 the spectrum is pinched with suppressed low and high energy tails, whereas for 
Pi < 3 the spectrum is anti-pinched (broader). The neutrinos decouple at a smaller radius 
and will therefore be hotter. Due to different charge-current interaction also the electron 
antineutrinos will decouple before the electron neutrinos. Correspondingly, a hierarchy of 
the form {E^) < {E^} < {E^) is expected. 

The unoscillated flux at distance D from the SN is given by 

n-^.^nE), (2) 

and we will flx the distance to 10 kpc. The neutrinos which are emitted from the interior of 
the star may undergo various flavor transitions due to MSW matter effects, when passing 
through the outer layers of the supernova. For the normal hierarchy, the the z/g and Ue 
survival probabilities. Pee and Pgg, respectively, are approximated by 

Pee ^ PH\Ue2\' + {I - PuWesl^ (3) 
P,, ^ \Uel\' , (4) 
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where Ua,i ,a = e,fi,T ,i = 1,2,3 are the Maki-Nakagawa-Sakata neutrino mixing matrix 
elements and we have used here the standard parameterization. We have fixed the value of 
the solar mixing angle 9i2 = 0.575 rad (sin^ 9i2 = 0.3) and the atmospheric mass squared 
difference Amg^ is set to be 3 x 10"^ eV^, as their variation have little impact on the final 
neutrino fluxes. We will disregard Earth matter effects and therefore the exact value of solar 
mass squared difference is irrelevant. Also the value of the atmospheric mixing angle is not 
relevant as the muon and tau neutrinos are indistinguishable. Ph is the hopping probability 
that can be written as 



Ph = exp 



1.08- lO^A'/' /|Am2 



sm ^13' 



E J \ 10-3 
The final fluxes arriving at Earth are simply given by 



41/3 



(5) 



F,^ = F°Pee + F°Jl-Pee), (6) 

Fp, = F9p,, + f9(1-P^^)^ (7) 



F.„ + F.. = F," fl - P..) + F^"^ (1 + 



Fp^ + Fp^ = po (1 - P,,) + Fl (1 + Pge-) . (9) 
These neutrino fluxes depend on our 10 dimensional parameters space, which is given by 



(Ei) , Li , Pi , i = e,e,x (10) 

13 



W^l, (11 



and we refer to the first 9 parameters as the SN parameters. 

There are four known detection channels for a water- Cherenkov detector: the dominant 
inverse beta decay (IB); the charge current on oxygen (CC-0); the neutral current on oxygen 



(NC-0) and the elastic scattering on electrons (ELAS). For a detai 



four possible channels in a Water- Cherenkov detector, please see 15|] 



ed description of these 



We will analyze three different cases for fitting the supernova and neutrino parameters. 

• Case A: We assume that only the inverse beta decay channel is available. In this case 
there is no dependence on the parameters. Moreover, since we consider only the 
normal hierarchy the allowed parameter space will even be independent of the Chooz 



angle (^13), since Pgg is constant (see eg. Fig.l of Ref.[15|])- Therefore, in this case 
only 6 parameters are left to be determined by data. 

Case B: Here we will assume that there are four available detection channels and the 
detected neutrino fiuxes are sensitive to all 10 parameters. 
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• Case C: In this case we also assume that all four detection channels are available, 
but we impose the constraint Lg = Lg, leaving us with 9 free parameters. 

The enormous flux of neutrinos which will arrived at the Earth from a galactic supernova, 
makes case A a highly pessimistic scenario, as it seems very likely that it will be possible 
to measure and also separate at least some of the other detection channels. Our main 
motivation for including this case is in order to be able to compare to previous works. 
Similarly, the condition Le = Lg has been used in previous works, so we include this case for 
easy comparison. We would like to note that we assume that all the four channels can be 
completely separated, like it was done in Ref . [3] • This is certainly not a realistic assumption 
but it allows us to establish what would be the best attainable results for a future Water- 



Cherenkov detector. We should however point out that with the addition of gadolinium [23 1 
it should be at least possible to do a fairly good separation by using the known directional 
forms of each event type. 

We simulate the signals for a given set of input values of the 10 parameters of Eq. f fTTl) 
and try to find the allowed regions of these parameters minimizing a function. For this 
purpose we define the simple function 

^bin ^jyth _ jY0bs-j2 
1 = 1 ' 

where N°^^ and N^^ are, respectively, the simulated and fitted number of events in the i-th. 
energy bin, and we take the number of bins, A'^biru to be 40. All bins have a width of 2.5 
MeV and we set the threshold at 5 MeV. For case A, the function only includes events 
from the inverse beta decay channel; 

Xcase A = XlB i 

(13) 

whereas we use all four available channels for cases B and C 

XCaso B, Case C = XiB + Xcc-o + Xnc-o + Xelas ) (14) 

and we calculate the 3 a allowed areas in the projected two dimensional space using 2 degrees 
of freedom. 

Naturally the allowed areas will depend on the assumed true values of the parameters, 
what we will refer to as the input parameters. We will just look at one illustrative example 
and for that we have chosen the input values (Ef.) = 12 MeV , (Eg) = 15 MeV , (E^) = 
18 MeV and furthermore we take all the integrated luminosities to be equal; Li = 0.5 x 10^^ 
Ergs. The spectral indices are assumed to be /3e = 5, /3g = 5 and Px = 4. Finally, the 
neutrino parameter sin^ ^13 = 10~^. These values are all chosen equal to the input values of 
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Ref. 171] as this allows for an easier comparison. For the detector size we have also taken the 



same value as in Ref. [17| with a fiducial volume of 720 kton. This value is somewhat high 
but it should be remembered that the fiducial volume for a galactic supernova detection 
is larger than that for eg. atmospheric neutrino detection and therefore this value is not 
unrealistic. 

based code for finding the Sex allowed regions. A number of 



We will use a MINUIT 



constraints on the parameters have also been implemented; 

5 < (£'e)/MeV < 17, (15) 
5 < {Ei)/MeY < 35 , i = e,x (16) 
< A < 25, , i = e,e,x (17) 

such as to avoid extremely unrealistic values. Also, for the case A, the ratio between the 
neutrino fluxes are not allowed to become too large. Especially this means that we take 
the minimum value to be the minimum which is the local minimum closest to the input 
values. To be specific, we set all parameters to their input values and let MINUIT fall into 
the local minimum from there. 

III. DISCUSSION 

In this section we will present our results, and also compare them to previous works. Fur- 
thermore, we investigate ways to break down the degeneracies which exists in the parameter 
space. 

In FigH] we show the 3 a allowed areas for the three cases A, B and C in six different 
planes of two parameter projections. It is seen that the allowed regions are quite different 
in each case. The allowed regions for case A are the largest and the ones for case C are the 
smallest, which is clear from the definitions of the three cases. The only exception is for /3e 
which remains completely undetermined even in case C. 

By looking at the contours for case A, we see that we overall agree with the findings of 
Ref. 171]. Indeed the degeneracy in the {E^) - (3x - parameters is also clearly exhibited in 
FigH] (panels 3 and 5). Although, there are minor differences in the size of the areas, the 



shapes of the allowed areas are the same as in Ref.[17j. The small differences must be due to 
slightly different analysis procedures. For instance, the width of the energy bins are different 
in the two analysis. As already explained there is no sensitivity to the parameters as well 
as to sin^ ^13 in case A. The determination of the Ue parameters are not impressive but still 
fairly good. The uncertainty on the value of {E^) is about 5 %, Lg is determined within 50% 
and /3e is consistent with values in two distinct island, one including the input value. Clearly 
the most important and problematic fact is that the parameters are left undetermined. 
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FIG. 1: The 3 a allowed regions for the three cases A, B and C as explained in section [III The 
outermost dark (blue) contour assumes that only the inverse beta decay channel is available, and 
thus only depends on 6 variables (case A). The middle dark (red) contour assumes that all 4 
channels are available and in this case there is in total 10 free variables (case B). For the light 
(cyan) contour we also assume that all 4 detection channels are available, but we implement the 
condition Lg = Lg, thus leaving only 9 free variables (case C). For the case A, there is no dependence 
on the electron-neutrino parameters and therefore this case does not constrain any of the variables 
in frame 1. 
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The degeneracy extends over the whole reahstic range of parameter. Therefore, without 
further input, it is unfortunately impossible to say something conclusive about the initial 
flux. 

In case B we find rather large allowed regions for some of the parameters, but considerably 
smaller than for case A. Although the degeneracy is partially broken, there is still correlation 
among the Vx parameters so a large elongated area is still seen in panels 3 and 5 of Fig{TJ We 
believe that the main reason for the fact that much of the degenerate area survives is simply 
that the other detection channels (CC-0, NC-0 and ELAS) have a much smaller number of 
events. For an idea of the orders of magnitudes of each event type, please see Table II in 



Ref. [l5l]. By comparing with the results of Ref.[15|] it becomes clear that the inclusion of 
the pinching parameters have a very large impact on the size of the allowed areas for the 
average energies, in particular that of the Vx neutrinos. This is also illustrated in Fig|2l The 
fact that, to a large extend the degeneracy is still left even after including all four detection 
channels, should be viewed as a large obstacle for the determination of the Vx parameters. 
Below we will study some possible ways of overcoming this degeneracy. 

For case B it is seen that the determination of the value of sin^ ^13 is quite good. It should 
be remembered that for values of sin^ 6*13 < 10"^ there is no dependence on this variable as 
the H resonance is always non-adiabatic. Furthermore, the values of the Ve parameters are 
somewhat better determined than for the case A. The value of {Eg) is determined within 
about 3 % and /3g is constrained within about 8 %. The electron-neutrino parameters are 
rather difficult to determine as there is not enough statistic in the detection channels. 
Nevertheless, in view of this, the allowed area of {Ef.) is rather small, {E^) is determined 
within 13% but one can give only a lower limit for the luminosity L^. > 0.1 x 10"^^ ergs. 

Let us next look at the contours for the case C in Fig{Tl These contours are in general 
very small. In particular, the degeneracy in the Ux variables is broken down and the averaged 
i>x energy is determined within 5 %, which is a big improvement when compared to case B. 
Also the constraints on (Eg) is somewhat improved, whereas the determination of 613 and 
/?e have not really improved as compared to the case B. In contrast with the other cases, all 
the luminosities are constrained within about 10% for the case C. 

Of course the assumption of identical electron neutrino and electron antineutrino in- 
tegrated luminosities can be criticized as it is not based on a strictly physical condition. 
Although the two luminosities are expected to be at least of the same order of magnitude. 



the ratio of the two is still varying a lot in different SN simulations. In Ref. ISj this assump- 



tion was made ^ in order to decrease the number of free parameters. Even if this assumption 



''' In Ref. 15 1 also the pinching parameters were fixed and the fiducial volume was 75% smaller, i.e. smaller 



by a factor 540/720. 
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seems rather innocent we see that it has a large impact specially on the allowed region for 
the fx parameters. In fact, from Figl2] it is clear, that either the assumption, Lg = 
or the fixing of the /5j's, significantly diminish the allowed average neutrino energies, most 
pronounced for {E.j.). One can also observe, by looking at panels 3 and 5 of Fig. [H that 
the fixing of I3x alone would allow one to determine quite well all the mean energies and 
luminosities (as long as the total luminosity is independently constrained). 




15 17 18 19 20 21 22 / ' 1 ,q-7 , p-^ ^q-i -^r^-^' 

<E^> (MeV) sin ^ 6 



FIG. 2: The innermost dark (black) contour shows the 2>a allowed region when fixing the /Jj's at 
their input values. We assume that all four channels are available. The outermost dark (red) and 
the light (cyan) contours are identical to the ones in Fig{T]and shown for easy comparison. 

However, neither the constraint Lg = Lg nor the fixing of the pinching parameters, 
influence the allowed region for the Chooz angle in a sizable way. In the left panel of Figj2l 
we have explicitly shown this by calculating the contour obtained when fixing the /3j's at 
their input values. Therefore, a nice feature that is kept even when including the pinching 
parameters is the rather fine determination of the Chooz angle as seen in the last panel 
of Fig{TJ This is particularly important, since for very small 6'i3, i.e. sin^ 6*13 H 10~^, the 
observation of a galactic supernova is, at least at present, the only way of determining (or 
constraining) 9i^. 

Another important question that arises is whether the hierarchy can be determined by 
the observation of a galactic supernova? Here we consider simply the observation of the 
neutrino fluxes without the measurement of Earth matter or shock wave effects. Let us 
briefly comment on this point. In this work we have fixed to neutrino mass hierarchy to 
normal. But, let us in the following discussion assume that the collective neutrino-neutrino 
effects a.e negligible even for the inverted nent.ino^ ma. Me.a.chy. In Ret. Q a discn^ion of 
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the prospects for determining the neutrino mass hierarchy is presented. In this work the /5j's 
are fixed and also the condition Lg = has been implemented. Despite these assumptions 
it seems clear to us that the results presented in Ref. 15| about the determination of the 
neutrino mass hierarchy will be valid even for a 11 dimensional parameter space (i.e. that 
of Eq. dlip along with the neutrino mass hierarchy). The hierarchy will influence the total 
number of events in each channel for large values of sin^ 6*13 (above 10~^) (see Table II in 
Ref. 1^ ) . Therefore the hierarchy can be determined by the combined use of all four channels 
in a Water- Cherenkov detector. The degeneracy that we have observed in the parameters 
on the other hand maintains the total number of events in each channel almost flxed and 
only changes the spectral form at low energies. Henceforth, we expect that for large values of 
sin^ ^13, the neutrinos mass hierarchy can be determined with very large statistic. Moreover, 
in the case of the spectral swapping due to collective effects, as described in Ref. [3], it 
might even be possible to determine the hierarchy for small values of sin^ ^13. 

One can wonder to what extend our results depend on the particular input we have 
chosen. Especially if the average energies are further apart, like eg. {E^) ^ {E^), would the 
degeneracy in the parameters be much milder? Such investigation are left for a future 
paper. 

In the following we will discuss ways to break down the (-Ex), (^x and degeneracy. As 
the main problem is pinning down the value of (-Ex), we will focus on the third frame in 
FiglH One possibility is that some robust features for the supernova parameters will emerge 
from supernova simulations in the future. These features can then be safely implemented 
in the fltting procedure and used to constrain the parameter space. This can be thought 
of as a standard supernova model, in analog to the use of the standard solar model when 
studying solar neutrinos. If a standard supernova model will be developed one can include a 
penalty in the function when eg. the ratio Lf./L^ differs from its central value as predicted 
by this standard model. We do not have such a standard model at our disposal yet, so to 
demonstrate our point we study some simple cases. 

Clearly, if the supernova simulations could tell us rather precisely the value of the pinching 
parameters, then this degeneracy would be broken. It should be remembered that we flt to 
the total number of events in the cooling phase. A slightly time varying spectral index, will 
cause the time integrated spectra to be broader. Therefore, we can expect lower values for 
the /3j's in the time-integrated spectra. By combining the information of the panels 3 and 
5 in FigHlone can directly relate constrains on (3x to constraints on {Ex). Notice that if we 
take the range value of 3-6, which is presently suggested by the SN simulations, then the 
constraints on the average energy is not improved for the case B. 

An appealing possibility is to use ratios of the individual luminosities. Let us deflne two 
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FIG. 3: The 3cr allowed regions for various assumptions about constraints on the supernova lumi- 
nosity ratios. We assume that all four channels are available in all the contours. The outermost 
dark (red) and the second smallest (cyan) contours are identical to the respectively case B and 
case C, in Fig(T]and shown for easy comparison. For the second largest (purple) contour we have 
constrained 0.5 < < 2.0 and for the third largest (green) contour 2/3 <Ci^3/2, where i = 1,2. 
The smallest (black) contour is calculated assuming ^2 = 1- 



One reason that these can be interesting parameters to use, is that there is a good chance 
that these can be constrained (in a model independent way) from supernova simulations. 
The physical processes that are behind the production of the neutrino flux are related to 
known physics. Nevertheless, there are a number of processes, such as nucleon-nucleon 
bremsstrahlung (A^A^ NNuu), neutrino- antineutrino annihilation (eg. z/gi^e vv) and 
various scattering reactions between neutrino and antineutrinos of different flavors, which 
of course complicates the calculation of the neutrino emission. From the present supernova 
simulations it seems that these ratios can at most be two and should be larger than one- 
half. It is important to notice that when increasing the value of {E^) within the degeneracy 
area, the ratios as well as ^2 also increases. In general, for values {E^) < 18 MeV the 
two ratios will be below one and for (Ex) > 18 MeV the ratios will be larger than one. In 
Figl3] we illustrate how the constrains on these ratios affect the allowed region. The second 
largest area in Fig|3] corresponds to the suggestions from present supernova simulations and 



ratios 



Lx 
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a smaller allowed area is obtained. In FiglS] we also show the contours when fixing ,^1 = 1 
(the same as the case C of Fig{T]) as well as ^2 = 1- In these cases a very good determination 
of {E^) is possible. In conclusion we find that constraining the luminosity ratios can be 
helpful for pinning down the supernova parameters. 




16 17 18 19 20 21 22 23 2-1 
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FIG. 4: The dashed dark (black) contour shows the 3o" allowed region when fixing the total 
luminosity of the supernova to the input value of 3 x 10^'^ ergs and leaving all parameters free. 
We assume that all four channels are available. The outermost dark (red) and the light (cyan) 
contours are identical to the ones in Fig(T] and shown for easy comparison. 



Next, we will look at the possibility that the total energy liberated in the supernova, which 
we will refer to as L^^ti has been already determined by some other method. In Ref. 22| it is 
suggested that a measurement of the electron neutrinos from the neutronization burst can 
be used to determine the distance to the supernova. This along with a optical observation 
of the SN, might be used to predict the total energy liberated. In FigH] we present the 
results for fixing Ltot to its input value. This reveals that even in this extreme situation 
the improvement in the allowed region is minor. Indeed, the neutral currents channels, the 
NC-0 and in parts also the ELAS channel, are already constraining the total luminosity. In 
case B the total luminosity is determined to about 4% accuracy (assuming that the distance 
is known). Clearly constraining L^^t is not very helpful for determining the supernova 
parameters. In fact, one might even expect that the fitting procedure (case B) will provide 
the best determination of Ltot • 
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IV. CONCLUSION 



We have shown that the inclusion of all supernova and neutrinos parameters is important 
for determining the allowed regions that can be obtained from the observation of the neutrino 
burst from a galactic supernova. We use a total of 10 parameters ((-Ee), (-Eg), {E^), L^, Lg L^, 
I3e, jSg, (3x and sin^ 9i^) and four detection channels (IB, CC-0, NC-0 and ELAS) that can 
be observed by a Water- Cherenkov detector to fit to such an observation. The degeneracy 



between (-E^), f^x and when using only the inverse beta decay channel [l7j, is only 
broken mildly by the inclusion of the other channels in a Water-Cherenkov detector. This is 
mainly due to the fact the number of events from these other channels are at least one order 
of magnitude smaller, and thus, in principle, the degeneracy could be broken by including 
more channels if statistics was not a limitation. Unfortunately, the supernova parameters are 
very difficult to determine due to this degeneracy. Especially, the parameters cannot be 
properly identified. We have discussed ways that supernova simulations can help overcoming 
this problem. A particular good way, seems to be to constrain the ratio of the integrated 
luminosities of the neutrino fiavors. 

We have demonstrated that the so-called Chooz angle, ^13, can in principle be determined 
very well even when freely varying all parameters (including the pinching parameters). In 
fact, whether or not the pinching parameters are freely varied, does not infiuence much 
the allowed ^13 region. This is so because the Chooz angle infiuence the ratios of the total 
number of events in each of the four different channels, whereas the degeneracy caused by 
the pinching parameters maintains the total number of events almost intact. In the same 
way we expect that the neutrino mass hierarchy can be determined for large value of sin^ ^13 
by a galactic supernova. 

On the other hand, data can only constrain the SN parameters describing the Vg. flux (/5g, 
{Eg) and Lg) and {E^. Without extra assumptions on the luminosities (ratios and/or total) 
or fix one cannot determine Lg, Lx and {Ex), even if the four detection channels could be 
completely separated. 

To summarize, the neutrino parameters can be determined quite precisely, whereas it is 
more difficult to determine the supernova parameters. 
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